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ABSTRACT 



A new formula is derived that generalises an earlier result for the infinite integral 
over three spherical Bessel functions. The analytical result involves a finite sum over 
associated Legendre functions, Pj^ix) of degree I and order m. The sum allows for 
values of I m| that are greater than I. A generalisation for the associated Legendre 
functions to allow for any rational m for a specific I is also shown. 
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1. Introduction 

The calculations of infinite integrals that involve a product of spherical Bessel 
functions have been the focus of many papers, amongst which are references [1-22]. In 
particular, reference [12] showed a derivation for an analytical evaluation of an infinite 
integral over three spherical Bessel functions, given by the form 

oo 

HX.M,Xr.k.k.k.)^Jr'nA^MM.rUr. (1.1) 



This type of integral has received considerable attention due to its use in nuclear 
scattering theory [9-14]. In this paper, we extend this result to analytically evaluate 
infinite integrals of the form 

00 

/(A; Ai, A2, A3; /ci/c2A;3) = J r^^^ jxAhr) j\2{k2r) jx^+xiksr) dr. (1.2) 



This type of integral has in the past been attempted both analytically [8] and 
numerically [17]. The authors believe that the evaluation carried out here is more 
compact, easier to work with and extendible. The new result, like the earlier one when 
A = 0, only works when |Ai — A2I < A3 < Ai -|- A2, i.e. when the integer indices Ai, A2 
and A3 satisfy the triangular condition, and Ai -|- A2 -h A3 is even. However, unlike the 
earlier result where the triangular condition forces ki, k2 and k-^ to form the sides of a 
triangle, this generalised integral can have non-zero values when the /c's do not form the 
sides of a triangle. Our new result is only derived for the case the A;'s do form the sides 
of a triangle. We also show the general result, i.e. for any values of k\, k2 and ^3, when 
A = 1. This new result involves finite sums over the associated Legendre function, 
P[^{x), of degree I and order m, whereas the earlier result involved sums over the 
Legendre polynomial, Pi{x). The sums can involve values for \m\ that are larger than 
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I. An extension for the associated Legendre functions which produces a formula that 
allows for such values is shown in appendix A. An application is shown to evaluate an 
integral involving four spherical Bessel functions and is used in the angular integration 
of the homogeneous and velocity isotropic Boltzmann equation. This integral, which is 
highly oscillatory and needs special treatment when evaluated numerically, is reduced to 
an integral over an associated Legendre function and a Legendre polynomial combined 
with algebraic factors which can easily be evaluated numerically. A comparison is shown 
in appendix B between our analytical result and the analytical evaluation carried out 
in reference [8]. 



2. Generalising the Integral Over Three Spherical Bessel Functions 

An earlier result [12] showed that an infinite integral over three spherical Bessel 
functions can be written as 



X 



M c l^ M M A3 

' oj \ c M-c I ^ ^ ' 

where A = {kj + - /c|)/2/ci/c2 and /3(A) = 0(1 - A)0{1 + A) with 9 the Heaviside 

function in half-maximum convention. Pi{x) is a Legendre polynomial of order I, 

Ai A2 A3 \ I Ai A2 A3 I 

is a 3i symbol and < > is a 6i symbol which can be 

OOOy \£A3-£/j 

found in any standard angular momentum text [23, 24]. Note that the summand in I 
vanishes unless |Ai — (A3 — C)\ < I < Ai + A3 — £. Now, multiply both sides by 
Ajg^"*"^ and integrate over from to where K can have any positive value. The 



left hand side involves the integral (see [25], eq. 5.52-1, page 661) 



K 



J kt^^jx,{k3r)dks = ^K^'+^jx,+i{Kr). (2.2) 



The right hand side involves the integral 



K 



J = j piA)k3PiiA)dk3. (2.3) 







Now, ksdks = —kik2dA. So, when A;3 = then ki = k2 and A = 1. Also, when 
ks = K,A = A', where A' = {kj + kj - K^)/2kik2. Hence 

1 

J = kik2 J PiA)Pi(A)dA. (2.4) 

A' 

If A' > 1, i.e. K <\k2- ki\, then J = 0. If -1 < A' < 1, then 

J = A;iA;2^(A0(l-A'')V2p-i(A'), (2.5) 

using [15] 

1 1 

p-^(a;) = (1 - x2)-"^/2 J ...J Pi(x) {dxr, (2.6) 

X X 

where PJ^ix) is the associated Legendre function of the first kind of degree I and order 
m. If A' < -1, then 

1 

J = kik2 e[K - {ki + k2)] J Pi{A) dA = 2kik2 e[K - {ki + k2)] 0- (2.7) 



-1 



Hence, the result is 



^ f Xi Xs-jC l\ f X2 jC l\ ( Xi X2 A3I 

Y vo oy\oooy\£A3->c/j' 

+ (-1) 2 (2Ai + 1) (2A, + 1) [2X2) ^["^ - + ^'^^ 

(2.8) 

The result (2.8) applies for any values of ki, k2 and K. The second term on the 
right hand side of (2.8) is also consistent with equation 6.578-4 of reference [25] when 
A3 = Ai + A2. This equation applies for K > ki + k2 and can be written in the form 

00 

(2.9) 

k^' k^^ r(A3 + 3/2) 

^ r(A2 + 3/2) r(Ai + 3/2) ■ 

If we assume that K satisfies the triangular condition, then the second term in (2.8) 

vanishes, and by repeated integration one arrives at (renaming K back to ^3) 

/(A ; A. A. A3 ; *3) = ' (2A3 . 1)'/^ (I) " 



X 



X 



(2.10) 

where A, Ai, A2 and A3 are all greater than or equal to to ensure convergence of 



the integral. Equation (2.10) is a new result that generalizes the infinite integral over 
three spherical Bessel functions. In general A can exceed the value for I. In this case 
the definition of Pp{x) needs to be extended for \ m\ > Z in a way compatible with 
eq. (2.6). Appendix A discusses such a generalization for the associated Legendre 
functions. Several explicit results for eq. (2.10) are presented in appendix B. For some 
small values of the indices these have been compared with a different existing result for 
the integral as will also be discussed there. 

3. Special Cases and Identities 

Equation (2.10) can be reduced to a known integral [3] for the case A3 = 0, i.e. 
Ai = A2 = A' with the result 

/(A; A', A', 0; k.kM - 3^ d - AY/^P-(A). (3.1, 

when ki, /c2 and ks satisfy the triangular condition 

\ki - k2 \ < k3 < ki + k2. (3.2) 

If we set A = A2 = and Ai = A3 = A' in eq. (2.10) and equate it to eq. (3.1) after 
setting A = and interchanging /c2 and k^, the following sum rule for the Legendre 
polynomial is obtained 

m E Q) (- 1)' PcM = m {ff (3.3) 

where ry = (fc^ + fcf - kl)/2kik2 and rj' = {kf + k^ - k2)/2kik3. Equation (3.3) is 
consistent with the result of reference [26]. In a triangle of sides ki, k2 and k^, rj is the 
cosine of the angle facing side ks and 77' is the cosine of the angle facing side k2. 
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4. Applications 



One application is an integral over four spherical Bessel functions which arises in 
the angular integration of the homogeneous and velocity isotropic Boltzmann equation 
[19,20] 

I(L; N + M - L,0, iV, M; h k2 ks k^) 

7 _ (4-1) 
= j "^^ ^ jN+M-Likir)jo{k2r)jN{k3r)jM{k4r)dr, 





where L, N and M are non- negative integers with N + M > L and ki, k2, ks and ki 
form the sides of a quadrilateral. Using the closure relation for the Bessel functions 
[2,19], this integral can be written as 



^ J K'^dKI{L; N + M -L, 0, N + M - L; kik2K) 



X /(O ; TV, M, TV + M ; /C3 /c4 K). 



(4.2) 



Now, using eq. (2.10) 



N+M-L 



I(L;N + M-L,0,N^M-L;k,k,K) = ^ (|) 



(4.3) 



where ^ = (kf + k'2 — K'^)/2kik2 and using 



2/y/2 n i-v i'\ _ {-1)1 n 



(4.4) 
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The resulting integral is 



I{L; N + M - L,0, N, M; ki k2 ks k^) = 



lQkik2k'^k4^ 

I l\l M l\l -L \ 

X {ki k3f+^ k^^ ^/2{N + M) + l 



N M N + M 


N+M-LN+M _ /OAr^OA#\ 1/2 



^ ,2N + 2MY^' fN + M-L\ 

£=0 C'=0 ^ ^ ^ ^ « 



X 



N N + M - C l\ / M C l\ ( N M N + M 

ojlo 0)\C'N + M-C' I 



X S{kik2k3k4;LMNC'l), 

(4.5) 



where, using ^' = (/cg + /c^ - K'^)/2k3k4, 



oo 

^(/Ci/C2/C3fc4; LMiV£'/) = 1 (1 - C')"^/' i^r^'^lO ^/(O rfi^, (4.6) 



is an integral that is easily done numerically as it converges rapidly. 



5. Conclusions 

A new generalised analytical formula for the infinite integral over three spherical 
Bessel functions was shown. It involves finite sums over the associated Legendre 
function combined with angular momentum coupling coefficients 3j and 6j symbols. 
The sums involved values of the order which exceeded the degree, which is in confiict 
with the usual definition of the associated Legendre functions. An extension for the 
associated Legendre functions is shown in appendix A. 
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APPENDIX A: A Generalised Solution 
to the Associated Legendre Equation 



The associated Legendre function of the first kind, PJ^ix) is a solution of the 
differential equation [25] 



(1 _ - 2x^^ + m + 1) - Prix) = 0, (A.i) 



given by 



Pr{x) = (-1)^ ^ •;,/ -^^j— (a;2 - 1)', M.2) 



defined for integer degree I > and integer order m (with |m| < Z). In this appendix, it 
will be shown that another solution exists for this differential equation that extends the 
associated Legendre function to any rational \m\ (|m| < oo) for integer I. This solution 
is irregular at a; = 1 when m > I and irregular at a; = — 1 for m < I. The formula can 
be used to derive closed form expressions for P^ for a specific / and any rational m. 

If we set Z = in eq. (^-1), it is easy to show that the following function, Pq^{x), 
is a solution: 

/ 1 _ 7- \ — m/2 

P^ix) = (—^) , (^.3) 



1 + x 

where coefficient that depends on m. Now, if we assume that the dependence 

on I can be separated, then in general Pf^ can be written as 

Pr{x)^AUx)P^'^{x), (A4) 

where Ai^j^ix) is a coefficient that depend on I, m and x. To determine the coefficient 
(up to a constant), substitute eq. (^.4) back into eq. (^.1). The resulting differential 
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equation for A; „j is 

(1 - x') ^-^^ + (2m - 2x) + 1(1 + l)AUx) = 0. (^.5) 

This identifies Ai^^n ^ the Jacobi polynomial, [25] -P^~"^' "^{x), defined by 

Hence, equation (A. 4) becomes 

■1 -a;\-W2 



where bi^m is a constant that depends on I and m. Using the normalisation for the 
associated Legendre functions 

-1 

and the normalisation for the Jacobi polynomials 
-1 

one can square eq. {A.7) and integrate it over x from —1 to 1, then compare it with 
eq. {A.8) to obtain an expression for bi^^: 

\l — m\l 

where the absolute value in |Z — m|! is introduced to allow m to go from — oo to oo, and 
the expression reduces to {I — m)\ for —l<m<l. Equation (^.7) then becomes 

which agrees with eq. (2.4), page 12 of reference [27] and is valid for integer degree 
I > and rational order — oo < m < oo. 
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When —l<m<l, equations (^-11) and {A.2) are equal, leading to the identity 

for —l<m<l. 

Another property that can easily be found is 



Pl-i-x) = -(x), (A.13) 



for positive integer Z and any rational m. 

Using equations (^.7) and (A. 10), one can express the associated Legendre function 
in terms of the Jacobi polynomial as 

Using the recurrence relation for the Jacobi polynomials [25] 

2(n + 1) (n + a + /3 + l)(2n + a + /3) P^+f ^(x) = (2n + a + /3 + 1) 
X [{2n + a + p){2n + a + p + 2)x ^ o? - ^^]Pt''^\x) - 2(n + a)(n + ^) {A.15) 
X {2n + a + /3 + 2)P^''_:f \x), 

we find the following recurrence relation for the generalised associated Legendre 
functions 



(l^ - w?) 



Using this recurrence relation and eq. (^-11) one finds the following closed-form 
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expressions for the associated Legendre functions at specific Vs and any rational m 

1 /I -1- rn/2 

(r^) • 

P^"(a;) = \2-m\ \ (3^^^ " 32;m - 1 + m^) (j^^j , (A19) 

P^{x) = — - {15x^ -15x^m-9x + Qxm'^ + 4:m-m^) (^^^^]"'^^, {A.20) 
3 — m ! \1 — xJ 



Pr(x) = — - (lObx^ -105x^171 -90 x'^ + 45x^171"^ + 55 xm 

I 4 - m I ! 

- lOxrW^ + 9 - 10 + m^^ 



(A21) 



APPENDIX B: Explicit Expressions 

Here we provide some explicit expressions for the integral (1.2) and compare to 
those obtained from a result given in [8] , which reads in our notation 

/(A; Ai, A2, Ag; /ci k2h) = - tt 

X (^sgnih + ks + k2)F^i^l^l{ku k2, ks) 
+ {-lf^sgn{-ki + k2 + /^3)i^A~j;(-^i' ^2, ks) (B.l) 
+ i-lf 'sgn{k, -k2 + k3)Fp,lll^-^{ki, -fcs, ks) 
+ {-if^+'^^sgni-k^ -k2 + /^3)Fa™J3 (-^1' -^2, ks)). 
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with 



m2 
m2+l 



pm,m,ms^, , > ^ ^ (Ai + mQ! (-1^ ^ (Aa + ms)! (-1)" 

''io(^3 + A-m3)!m3! (2/c3p+^ (mi + ma + m3 + A)! 

(S.2) 

Restricting ourselves to values of ki, k2 and k:^ which satisfy the triangular condition 
and using symmetries of Fx'x^x^li {.B.l) can be simplified to 

/(A; Ai, A2, A3;/ci, ^2, /C3) = -tt /3 (A) 

^ ^ (-Ip (Ai + mi)! ^ (-in (A2 + m2)! 

y^' (-ir(A3 + A + m3)! 1 

^o(A3 + A-m3)!m3! (2/c3p+M™i + ™2 + m3 + A)! ' ^ - ^ 

{k^ + A;2 + A;3)^i+"»2+m3+A _ ^_^^\.+m, ^^m,+m2+rm+\ 
_ ^_-|^-^A2+m2 ^^mi+m2+m3+A _ ^_-|^-^A3+A+m3 ^^mi+m2+m3+A 



with ci = {—ki + /c2 + ^3)) C2 = {ki — k2 + ks) and C3 = (/ci + /c2 — ^3)- For numerical 
computations this expression is more suitable than (-B.1) in which almost cancellations 
between the different terms including sgn functions can lead to large relative errors. 
For some small values of A, Ai, A2 and A3 this yields the following explicit results: 

/(0;0,0,0;/ci/c2/c3) = ^^^, 

4:kik2ks 

7(0; 1, 0, 1; kik2k3) = --^ ^!^^^\ ( - 2kik3 + ki^ + 2/cici - 2A;iC2 - 2A;iC3 + ks^ 

IQki k2ks V 

ey <~\ <~\ Q \ 

- 2A;3Ci - 2A;3C2 + 2^303 -k2 + ci + C2 +03), 



14 



/(O; 1, 1, 0; hkiks) = - ""f ( - 2kik2 + ki^ + 2/ciCi - 2kiC2 - 2kiC3 

16 kik2k3 V 

+ /C2^ - 2ci/c2 + 2A;2C2 - 2A;2C3 - ks^ + ci^ + C2^ + cs^) , 

/(I; 0, 0, 0; A;iA;2A;3) = 7^ ^^l fa^ - 2A;3Ci - 2/C3C2 + 2^303 - ki^ - 2kik2 

16 kik2k3^ V 

- A;2^ + Ci^ + 02^ + 03^) , 



/(I; 1, 0, 1; kik2k3) = 7^ ^2^/^^ (2^2 V + 4/c3'ci2 - Uki^^s + C2^ + 8A;iA;3'ci 

64 kik2kz V 

- 8A;iA;3^C3 - 8A;iA;3^C2 + 12A;iA;3C2^ - 12A;iA;3C3^ + &kik2^ - 4A;iC2^ - 4A;iC3^ + Akici 



+ 4A;3C3^ + 4^3^02^ + 4/c32c3^ - 4/c3C2^ - 4A;3Ci^ - 2ki^k3^ + 8A;i3/c2 + ci^ 
+ C3^ + 3A;i4-/c2^-A;3^), 

/(I; 1, 1, 0; kik2k3) = , \^}f} , ( - 4/e2Ci3 - 6k2V + 24kik2k3C2 - Ukici^s 

iyz ki k2 ks ^ 

+ C2^ + 12A;iA;3C22 - nkikac^^ - 12kici^k2 - I2k2k2,c2^ + 12A;2A;3Ci2 - 12A;2/c3C3^ 
+ 12kik2Cz' - l2kik2C2^ + 2Akik2kzcz + 2Akik2kzci - Ak2Cz - 6A;i^A;2^ - 4A;iC2^ 
- 4A;iC3^ + Akici^ + 4^303^ - 4A;3C2^ - Ak^ci^ - Qkik^ + 4A;2C2^ + ci^ + 03^ 
+ 3A;i^ + 3A;2^ + 3A;3^). 

{BA) 

Higher order expressions are not shown here for the sake of brevity. On the other hand 
we get from eq. (2.10): 



/(A; Ai, A2, A3; kik2ks) = 



2^+2(Ai + A2 + A3 + A)! (5.5) 
X kr^'^+'h2-^'^+'hs-^^'+'+'^AiX; Ai, A2, A3), 
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with 



^(0; 0,0,0) = 1 
^(0;1,0,1) = -2A;i(-A;i + A;2A), 
A(0; 1,1,0) = 2/C1/C2A, 
yl(l;0,0,0) = 2/ciA;2(l- A), 
^(1; 1, 0, 1) = -6 ki^k2 (-2 ki + k2A + k2) (1 - A) , 
^(1; 1, 1, 0) = 6 kiW (A + 1) (1 - A) , 
A(2;0,0,0) = 4 (A-1)2A;22A;i2, 
A{2- 1, 0, 1) = -16 (A - if (-3 /ci + 2 /C2 + /C2A) k2^ki^, 
A{2; 1, 1, 0) = 16 (A - 1)^ (2 + A) k2^ki^. 



To test our new result we have compared several explicit analytic expressions obtained 
from (2.10) and (-B.3) respectively (including the ones given above), and found full 
agreement. The advantage of equations (-B.5) over {BA) is their compact form in 
terms of generalized Legendre functions. This property will be exploited in [28] where 
we study the integral 



00 




(S.6) 
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